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A new well-posed nonlocal Perona-Malik 
type equation 

Patrick Guidotti 



Abstract 



A modification of the Perona-Malik equation is proposed for which the 
local nonlinear diffusion term is replaced by a nonlocal term of slightly 
reduced "strength" . The new equation is globally well-posed (in spaces of 
classical regularity) and possesses desirable properties from the perspective 
of image processing. It admits characteristic functions as (formally linearly 
"stable") stationary solutions and can therefore be reliably employed for 
dcnoising keeping blurring in check. Its numerical implementation is stable, 
enhances some of the features of Perona-Malik, and avoids problems known 
to affect the latter. 



1. Introduction 

The Perona-Malik equation [14] 

u t- di <TTWW Vu) = ° (L1) 

o 

on a domain J? c K™ is usually complemented with homogeneous Neumann 
boundary conditions. The nonlinear diffusion coefficient 

a(u) 



1 + |Vu| 2 

is sometimes replaced by a(u) — e~l Vll l 2 . Pcrona and Malik originally in- 



troduced ( 1.1 1 in the context of image processing with the aim of dcnoising 
a given image uq while at the same time controlling blurring. The latter 
is unavoidable if a linear diffusion equation is used instead. The nonlinear 
diffusion determined by a(u) significantly impedes diffusion in directions of 
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steep gradients which correspond to sharp edges in the image. If one consid- 
ers a large gradient across a surface, then the nonlinear diffusion coefficient 
leads to diffusion in tangential directions and backward diffusion in normal 
direction. This is the reason why models of this type are often referred to 
as anisotropic diffusion models. 



Numerical approximations of ( 1 . 1 1 have confirmed this prediction and have 



been successfully employed as a consequence. It has, however, been observed 



that discretizations of (1.1) often lead to stair-casing in the numerical so- 
lution (see Figure [T]), or, more in general, to the preservation of gradients 
which do not represent any image feature but might simply be due to the 
presence of noise (cf. [MlH] and Figure |1| . Equation ( |1.1| poses many a 
challenge from the analytical point of view. It is an example of so-called 
forward-backward diffusion. In |10j the author shows that it is not well- 
posed. To circumvent this problem a variety of regularization techniques 
have been proposed. Spatial regularizations were the first to appear in [5) 
The gain in analytical tractability is, however, offset by the introduction 
of blurring. Temporal and spatio-temporal regularizations have also been 
considered, see [T2"I8 7 . 4 . 5 . 2 J . In [5j the author considers a purely tempo- 
ral regularization which leads to time-delayed Perona-Malik equation which 
mimics the implicit linearization procedure common to all numerical dis- 



cretizations of (1.1 1, that is 



u 



n+l 



u n + div (a(u n )Wu n+1 )6t n>0, 
1 u° — uq n = , 

whereby the nonlinearity is always evaluated at the previous time step. The 
introduction of a time delay leads to a tractable analytical problem and 
does not seem to negatively impact the desired features of a corresponding 
numerical implementation. The author proves local existence of a regular 
solution for the delay equation. 

Variational techniques have also been utilized to gain analytical understand- 



ing of ( 1.1 ) in a one-dimensional context. The Perona-Malik equation can in 
fact be viewed as the gradient flow associated with the non-convex energy 
functional 

E(u)= I log(l + |Vm| 2 ) , ueH'ffl). 
Jo 

Using concepts and techniques related to Young-measure solutions various 
results have been obtained. [T5JCEHI prove instability of certain solutions and 
the existence of infinitely many solutions of a certain (weak) type. 



Equation ( 1.1 ) was motivated by its potential applications to image process- 
ing. Whereas it poses significant and interesting mathematical challenges, 
its form is clearly not dictated by any physical principle. A modification 
is therefore proposed in this paper which is globally well-posed, and allows 
for special natural functions to be stationary solutions. The latter leads 
to a desirable dynamical behavior from the practical point of view. It also 
leads to stable pseudo-spectral discretizations with satisfactory properties 
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from the perspective of image processing. It can also be viewed as a new 
regularization technique. It differs, however, substantially from other regu- 
larization techniques in that it actually makes piecewise constant function 
become stationary solutions and in that the degree of regularization can be 
more finely tuned. In particular it does not regularize at any given specific 
scale. 

To simplify the discussion, analytical considerations will be restricted to 
the one dimensional setting. A representative numerical experiment, how- 
ever, will be presented also in the natural two dimensional setting. Roughly 
speaking, the modification proposed consists in replacing the nonlinear term 
by 

a ^ u) = ; u \^-< 12 (1 - 2) 

1 + [(-A) — u] 2 

for s G (0, i) and where A is the Neumann Laplacian. The fractional powers 



appearing in ( 1.2 1 can be defined in various manners depending on the choice 
of function spaces in which the equation is considered. Since one is interested 
in denoising while preserving sharp edges, it is clear that a £ should provide 
the very same benefits for the kind of large gradients one is interested in 
preserving. The associated equation 

"*-(! + [(- V««]* U *)" = (L3) 

has, however, the advantage of being quasi-linear and of avoiding backward 
diffusion altogether. The analysis and the numerical experiments performed 
in this paper therefore also show that backward diffusion is not an essen- 
tial ingredient in order to preserve sharp edges. For analytical reasons the 



ideas just presented will be performed on an equivalent formulation of (1.1) 



leading to an equation of the flavor of ( 1.2 1 



2. The Problem 

In its one dimensional formulation, the Perona-Malik equation reads 

'u* ~ (t^tU x ) x =0 in (0,1) for t> 0, 
< u x = at x = 0, 1 for t > , (2.1) 

u = u at t = in [0, 1]. 

Introducing v(x) = Jqu(u) dy as a new unknown, it is easily seen that it 
satisfies 

'vt - j^rV xx = in (0, 1) for t > , 

< v(t, 0) = , v(t, 1) = /o uo(y) dy for t > , (2.2) 

M '-) = I u o(y)dy in [0,1], 
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since it is readily verified that the average of uq is preserved over time using 
(2.1 ). By further defining w(x) = v(x) — x JoUo(y) dy , x e [0, 1] , it follows 
that 



w t 



Y^-w xx = in (0, 1) for t > , 

w (t,0) = Q = w(t,l) fort>0, (2.3) 

w(0, •) = J u (y) dy - (•) f u (y) dy in [0, 1] . 



Conversely, any solution of ( |2.2[ ) leads to a solution of (2.1 1 by setting 

u{x) = v x (x), x € [0, 1] . 

The evolution equation is clearly satisfied whereas for the boundary condi- 
tions the following argument is needed. The function 



v(x) 



satisfies 



Vt 



u(y) dy, x e [0, 1] , 
u x (0) 

; v ™- i+^(o)- 



"x(0) 



Since v — v by v(0) — and since v satisfies (2.2 1, it follows that 1 ^_ tl 2 ( ) 
and therefore u x (0) = 0. Now 



implies that 



u(t,y) dy = v{t,l) - «(t,0) = 

_u x (l) 
1 



= / u t (y)dy 



Uo(y) dy 
u x (0) 







1 + «I(0) 

and therefore that it x (l) = 0, too. Notice that if u is a piecewise constant 
function, then v is a continuous piecewise affine function. Thus piecewise 
affine functions play the same role for (2.3 1 as piecewise constant functions 
do for (pEIL 



Equation (2.3 1 is fully nonlinear and clearly presents the same analytical dif- 
ficulties as the original Perona-Malik equation. The following modification 
is proposed 

ju t - a e {u)Au = in (0, 1) for t > , 
i u(0, •) = u in [0, 1] , 



(2.4) 



for 



a e (u) = (l + K-A^-V) ,ee (0,1/2), (2.5) 

and where A now denotes the Dirichlet Laplacian. This is perfectly analo- 
gous to ( 1.2 1. Notice that the intensity of the nonlinearity is barely reduced. 
It is therefore to be expected that solutions of (2.4 1 with e > behave sim- 
ilarly to the solutions of (2.4) with e = for which the original (1.1) is 
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recovered, at least at the discrete level. At the continuous level (2.1 1 is ill- 



posed and such a claim is not very meaningful. It is therefore likely that 



(2.4) can provide a viable practical dcnoising tool while preventing blurring. 



Remark 1. It is easily verified that 

dtj 2dt7 7o l+^-'u] 2 



for v — u x and a smooth solution u of (2.4 1 



3. Global Existence of Smooth Solutions 



It turns out that (2.4) can be shown to be well-posed in the classical 
sense as soon as e > 0. In order to prove this, some notation needs to be 
introduced. For a G [0, 1) let 

rafn u _|{«6C a ([0 ) l],K)Wi) = ) i = 0,l} a 6 (0,1), 

C ° (U ' >- |c([o,i],m) a = 0, 

where C a ([0, 1],RJ is for a > the standard space of Holder continuous 
functions endowed with norm 



The semi-norm [-] a is given by 

[u]a := sup u e C a ([0 

if a > 0. Given a S [0, 1), the operator A Q is given by 

fdomA a = C 2 +a (0,l) ■= {u e C 2+Q (0, 1) | u, u xx e C^(0,1)} for a > 
jdom A = Cq(0, 1) := {u G C 2 (0, 1) | = , j = 0, 1} for a = 

and 

A a u = , u G dom(A a ) . 
Once a G [0, 1) is fixed, the following simplified notation will be used 

A:=A a , E := C?(0, 1) , E t = C 2+Q (0, 1) . 

It follows from [11, Corollary 3.1.21, Corollary 3.1.32, Theorem 3.1.34] that 
A is a sectorial operator and therefore generates an analytic semi-group 
{e tA 1 1 > 0} on Eq. Then its fractional powers {—A) 1 , 7 G (0, 1), can be 
defined as the inverses of 

1 r°° 
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defined on their range, that is, by 

dom((-Ar) = R((-^)^)- 1 , (-A) 1 = ((-A)^)- 1 . 

The domain of (— A) 1 might not be an interpolation space between E\ and 
Eq but it satisfies the interpolation inequality 

\\(-Ayx\\ Eo < c\\xf-^\Ax\\l o , x £ dom(A) . (3.1) 

It can, however, be sandwiched between known interpolation spaces 

(E , £i) 7) i ^ D((-Ay) (Eq, Ei) Ji0O (3.2) 

where, for 7 £ (0,1) and p £ [l,oo], (Eq, Ei)~ iP is the standard real inter- 
polation functor. Observe that for a = it can proved that 



(£o,£i) 7 ,oo = C^(0,l), 7 G(0,l)\{-}. 



Holder spaces of the time variable will also be useful. For (3 £ (0, 1), let 

E :=C ([O^Eo)^ ■.= C 1+0 ([O,T],E o )n&([O,T],E 1 ). (3.3) 
If u £ Ei , then 

1 

;A 



\ + [{-Ay-tu{t)f 

is sectorial for every fixed t. By Holder maximal regularity (see |lll Propo- 
sition 6.1.3]) it follows that any solution of 

S V t~ l + [(-A)i-<u]* Av = ° fOTi >« 



V 



(0) = u Q 



satisfies u £ Ei for any uo £ E\ such that Auq £ (Eo,Ei)p !<x> . In order 
to apply the mentioned theorem, one needs also to take into account the 
additional facts that 

[* " ■= 1 + [( -A)i-e u (t)]^ J G &([0,T\,W,Eo)) 

and 

dom(A(i)) = dom(A(0)) , t £ [0, T] . 
It is now possible to prove the following result. 

Theorem 1. Let uq £ E± such that Auq £ (Eq, Ei)f} j00 . Then there exists 
T + (uq) and a unique u : [0,T + (uo)) — * E\ such that uL Q _ is a solution of 

(2.4) on the time interval [0,T] for any < T < T + (uq). Furthermore 
u£C 1+(i ([O,T+(u o )),E o )nC ([O,T+(uo)),E 1 ) 
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Proof. Let v € Ex and define <P(v) to be the solution of 

1 



;Au = , u(0) = Uq . 



1 + [(-A) 1 - £ v] 2 ' 
Define r = 2\\uo — <P(uq)\\-e 1 and let 

B r := {v eEi\ ||u-u |k < r} . 

It follows that 

PO) - u Q \\ Sl < \\<P(v) - #(u )|| El + II^K) ~ uolk = 

|| [d t - a s (v)Ay 1 (0 7 u o ) - [d t ~ a E (uo)A] 1 (0, u ) || Ei + ^ 
[fit - a £ ( Uo )A] _1 {l - [dt - a e (u )A] [d t ~ a £ (v )A) _1 }(0, u Q ) 



r 

Ex + 2 

< r 



if T > is chosen small enough since v € B r and i>(0) = uo, so that v is 
uniformly close to uq. To see this, consider 

\\[a £ (v) - a e (u )]^||£(E 1 ,E ) < IM W ) - a e("o)ll£(E ) 

<\\(- A y-%v-u )\\ mo (3.4) 



which follows from 



l + w 2 (t) l + w 2 (s) 



3 < IK S ) - w(t)\\ Eo 



by setting u; = (— A) 1 ~ £ v since w(0) = u Q and a e (uo) are independent of 
the time variable. To further estimate (3.4 1, one can use the interpolation 
inequality ( |3.1| and 

v(t) — v(s) — J v(t) dr and v S Eq (3-5) 

to obtain 

\\(-A)^[v(t) - u ]\\ Eo < c\\v(t) - uoll^lK*) - u ||| 

< crt^ 1 - 5 ^ 5 , s,t < T, 

and 

\\(-A)^[v(t) - v(s)]\\ Eo < c\\v(t) - v(s)\\ E -c\\v(t) v(s)\\ Eo 



< cr\t - sfi 1 -^ , s,t<T. 



It follows that 
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and finally that 

||[a e (t;) - a e (u )]^|U(E 1 ,Eo) < crT^ 1 -^ 

Thus the claimed inequality follows by a simple Neumann scries argument 
made possible by choosing T small enough. Next observe that — <P(v2) 

solves 

f [#(ui) - $(v 2 )]t = a e («i)A[^(ui) - <f>(v 2 )] + (a e («i) - o e (v 2 )) A<2>(v 2 ) 

\[*(«l)-*(«2)](0) =0 

It follows that 

- «P(«a) = [9 t - a e («i)^] _1 {0, [a e («i) - o e («a)] A<5(« 2 )} . 

Thus 

||*(«i) -^(« 2 )|k < c(r)||[a e ( Vl ) - a e (« 2 )]A^(w 2 )|| jEo 

Now, since i>i — w 2 S Ei it follows that 

\\A l - s [vx - u 2 ]|| Ed < c||«i - uallE.T^ 1 -^) 



by making use of the interpolation inequality (3.1 1 and of (3.5 1 in a perfectly 
similar way to the calculations used to obtain the self-map property. Thus 
Banach fixed-point Theorem implies the existence of a unique solution 

u € Ei 

provided the time interval is chosen short enough. It follows from [11 , Propo- 
sition 6.1.3] that 

u(t) e CE ,£i)/3,oo, * e (0,T] 

which gives Au{T) € {Eo,E\)f} t00 . Thus, by using the same argument, the 
solution can be extended to a larger interval of existence. Repeating this ar- 
gument indefinitely, a solution is obtained on a maximal interval of existence 
[0,T + (mo)) with the stated properties. 

Next classical Sobolev spaces are needed 

W£(0, l) = {ti£ L p (0, 1) | d j u G L p (Q, 1) for < j < 2} 

for pg (0,oo) where the derivatives have of course to be understood in the 
distributional sense. The simple embedding 

W;(Q > l)^C- 1 /f(0,l) > a P >l I (3.6) 

will be very useful in the proof of the next Lemma. 

Lemma 1. If it can be shown that u 6 L 00 ([0,T + (wo)),Wp(0, 1)) for all 
p € (l,oo), then the solution exists globally in time. 
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Proof. Observe that 

u - a e (u)Au € L oo ([0,T + (m )),L p (0, 1)) 

follows from the assumption and the form of a s . For s,t £ [0,T + (uo)] , one 
has 

\\A l -*(u(t)-u(s))\y 

<c\\A^(u(t) - u(s))\\ wlP < cllA 1 -^^) - u(s))\\ Lp 

< c \\A(u(t) - u(s)) \\i- e+p \\ fu{r) drf-f < c\t - " 



as follows from (3.6 1 and (3.1 1. Hereby it needs to be assumed that 

p = 2p - - > and < p < e . 
P 

This can always be achieved by choosing p large enough and yields 

A l -*u e C^"([0,T+(u o )],&(0,l)) . (3.7) 
Denote the function space in ( |3.7[ ) by Eo. Let then v be the solution of 

{' — a e (u)Av = , v(0) = uq 
on [0,T + (ito)]. It satisfies v S Ei and one has 

V \[0,T] =M l[o,T]' T<T+(u ) 

by uniqueness. Since v is uniformly (Holder) continuous, so must be u. It 
can therefore be extended as a solution to [0, T + (uq)\ and as a consequence 
of Theorem [l] beyond any finite time. The solution is therefore global. 

Lemma 2. Let uo S C°°([0,1]) satisfy compatibility conditions to all or- 
ders, that is, assume that uq S dom((— A) k ) , k S N. Then the solution of 
(2.4) satisfies 

ueC 1+ ^([o,r+K)),c~([oa])), 

for (} < e. 

Proof. Let u be the solution of ( |2.4[ ). Then 

(-Af-'u E C^[0, T+(u )) , C o a+2/3 (0, 1) 

since (3 < e and u is therefore also a solution of ( |1.1[ ) in Ei for a + (3. By 
repeating this boot-strapping argument indefinitely one obtains the claim. 

Remark 2. Time regularity can also be improved by similar boot-strapping 
arguments but it will play no role in this paper. 
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It is now possible to prove the following global regularity results. 



Theorem 2. The maximal solution of (2.4 1 exists globally. 

Proof. According to Lemma [l] it is enough to show 

\\Au{t)\\ W 2 <c,te [0,T+(u Q )) forpe (0, oo) . 

By Lemma [2] it follows from ( |1.1[ ) that 

j Aii - A(a £ (u)Au) = , 
I Au(0) = u . 



Thus v := Au satisfies 



v — A(a E (u)v) = , 
v(0) = u . 



(3.; 



Observe that a e (u(t))A- is, for any fixed time t, formally adjoint to the 
operator ^4[a e (u(i))-] . The evolution operator U^t^r) generated by A = 
a e (u)A on E satisfies 

|| U A (t,T)u h p (Q,l) < ||w ||l p (o,i) • 

This follows from the Trotter-Kato product formula (see [13]) and the fact 
that the operator 

c(x)A 

generates a contraction semigroup on L p (0,l) if ellipticity is assumed (cf. 
[5]). Let T < T + (u Q ) and consider the family of generators B = A(T - t), 
t £ [0, T]. Then an easy computation based on 

-^-J7 A (t,r)u = t)A(t)u , u G domA(O) , 

ur 

reveals that 

U£{T,T-t)v ,te[0,T], 



satisfies (3.8 1 if t/g on L p /(0, 1) is taken to be the evolution operator dual 
to Uk- Now, since 

l|C/*(^T)|| £(v) = ||C/(t,r)|| £(Lp) <l 

it follows that 

IK*)|| V < ||«o|| v = \\Au \\ Lp , 

for any p € (1, 00) and therefore 

\\Au(t)\\ Lp <c||u|| W 2 <c||Auo||l, 

for all p G (1, 00). The proof is thus complete since T < T + (uq) is arbitrary 
and the embedding constant c does not depend on T. 
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4. Stationary Solutions 



Among constant functions, u = is clearly a stationary solution of (2.4 1. 
It corresponds to a constant solution with value J Q uq(x) dx of the original 
Perona-Malik equation with initial value uq. It turns out, however, that 
other more desirable functions are equilibria of (|2.4[) as well. Even though 



the next result remains valid for (2.4 1, it is formulated in the context of the 



modified Perona-Malik equation with periodic boundary conditions instead 
of homogeneous Dirichlet conditions. This allows the simplify a little the 
argument. 

Proposition 1. Let u be a continuous l-periodic piecewise affine function. 
Then 

for A taken to be the realization of d xx in the space of Radon measures M w . 
The subscript n indicates that the periodic case is considered. 

Proof. If u has the properties stated in the proposition, it follows that 

n 

^xx^ ^ y Cj^Xj 
3=1 

for some constants Cj £ K and locations x\, . . . ,x n € [0, 1). The claim would 
therefore follow if it can proved that 

1 



is a continuous function vanishing at Xj , j = 1, . . . , n. It is clear that it is 
enough to consider the case of a function with a single kink. Furthermore, 
since periodicity is assumed, Poisson's summation formula implies that the 
singularity in (— A) 1 ~ e u in [0,1) is the same as the singularity at the same 
location for (— periodically extended to real line. By making use 
of standard localization techniques, it is therefore sufficient to consider a 
function of the real line with a single kink in the origin and A = d xx on the 
full line. Choose the prototype u(x) — \x\. Since 

(-A) 1 - e u=(-A)- e S 

it follows that JF((— yl) 1_e u) = and therefore that 

vl-s 



for C e — ^J~^T(l — 2e) sin(7re). Remember that it has been assumed initially 
that e € (0, 1/2) to make the nonlinearity strong enough. It follows that, in 
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the general case, (— A) 1 e u has integrable singularities at Xj, j = l,...,n, 
and is otherwise smooth by hypo-cllipticity. This shows that 

1 



1 + [{-Ay-^u} 



is a continuous function which vanishes exactly at xj, j — 1, . . . ,n. It is 
therefore a well-defined multiplier for c j^xj and one has 

1 

77, Au = 



as claimed. 

Remark 3. It should be observed that it is not possible (to the best of 
our knowledge) to make sense of ((OJ as generating a flow on M v (or on 
any of the standard Besov spaces containing Dirac distributions) and so 
the stationary solutions of the proposition are only formally such. A short 
computation would even show that they are formally linearly stable. This 



might help explain the fact that, in numerical calculations of (2.4 1, it is 
observed that solution starting close to such a stationary solution tend to 
stay in its vicinity for a long time before being driven to the trivial steady 
state. Even more is actually true. Numerical experiments (see Figure [2| 
show that smooth solutions tend to become piecewise affine at first but are 
eventually completely smoothed out. 



Remark 4. Since the weakened nonlinearity (2.5 1 does not give rise to 
forward-backward diffusion for the linearization, global well-posedness can 
be proved. Its linearization, however, will be closer and closer to the Perona- 
Malik linearization as e decreases. As a consequence more and more of its 
eigenvalues will become negative when linearizing in the presence of large 
gradients. This is supporting evidence that the numerical benefits of Perona- 



Malik should not go lost in the modified (2.4 1 



Remark 5. The new equation is globally well-posed for smooth initial con- 
ditions. It is, however, not immediately clear what their large time asymp- 
totic behavior should be. Numerical experiments seem to suggest that any 
solution will eventually converge to a trivial steady-state. This behavior 
might, however, be due to numerical diffusion. In the next section a reason 
will be given which seems to exclude the latter possibility. 

Remark 6. It is also observed that the original Perona-Malik equations do 
not admit piece-wise constant functions as formal stationary points unless 
an ad-hoc concept of generalized stationary solution is introduced (see [TO]) 
since the multiplication of distributions is in general not well-defined. Its dis- 
cretization, however, seems to exhibit "almost stability" of such solutions in 
the sense explained here. This might be an indication that the Perona-Malik 
equation is regularized by discretization and could explain the success of its 
implementations in spite of the ill-posedness of its analytical counterpart. 
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5. Numerical implementation and experiments 



Next a numerical discretization of (2.4 1 is derived in a periodic con- 
text and used to perform numerical experiments intended to illustrate and 
demonstrate the claimed improvements on the classical Perona-Malik equa- 
tion. The periodic Laplacian A is discretized spectrally by means of the fast 
Fourier transform F n 

A n = F~ 1 A n F n (5.1) 
where n = 2 m denotes the number of grid points used and 

D n = 4tt 2 diag((-| + l) 2 , (-| + 2) 2 , . . . , 0, 1, . . . , ^) . 

The time variable is discretized by forward semi-implicit Eulcr so that 

u k+1 = \id n ^ A n ] ~ V (5.2) 

where u k is the spatial n-vector at time k h t for the time step h t > 0. Observe 
that, setting e = 0, the classical Perona-Malik equation is recovered. 
The following experiments are presented here. First the evolution of the 
function 

u (x) := 100x 2 (l - x 2 ) , x € [0, 1] . 

is considered. Figure [l] depicts the first derivative of the function u k after 
100 time steps of size h t — 0.06 for m — 8. The blue curve corresponds to 
the Perona-Malik solutions and the stair-casing phenomenon is apparent. 
The red, magenta and cyan curves correspond to solution of the modified 
equation for e — 0.3,0.2,0.1, respectively. They clearly benefit from all 
effects of the Perona-Malik equation without leading to stair-casing. Some 
decrease in contrast is the only price paid. Notice that the black curve 
depicts the first derivative of the initial value. Comparison with the other 
curves reveals the de-blurring effect of the equation. It has been proved 
that continuous picccwisc affinc functions are steady-states for the modified 
Perona-Malik equation. There is, however, no mechanism by which a smooth 
solution can develop a singularity in finite time. The solution therefore feels 
the presence of these nontrivial steady-states but never leaves the regime 
where diffusion, slowly but surely, drives it to a trivial steady-state. This is 
exemplified in Figure [2] where various stages of the evolution of a smooth 
solution are with initial condition 

uo(x) = sin(27nr) + 2sin(47ra) , x G [0, 1] . 

is depicted. Clearly the convergence towards the trivial steady state might 
be due to numerical diffusion (see Remark [5]). However, Remark [l] provides 
an indirect way to test the numerics. If follows from Remark [T] that 

i rl rt pi 2 

u 2 x dx= / (u f x dx-2 / / f—^dxdr. 
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This conservation relation can be tracked for any smooth solution and, in 
the particular case considered, the relative numerical deviation observed 
between the left and right-hand-side of it amounts to a mere 0.2% for a 
fully converged solution. Scheme (5.2) also preserves the piecewise affine 
structure of initial values for large integration times. The solution with 
initial condition 



u {x) = 5- |10x-5| , x £ (0,1) , 
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Fig. 3. Behavior close to formally stationary solutions. 



is computed for 100 time steps with ht — 0.06 and m = 8. The solutions 
are plotted in Figure [3j The color coding is as in the previous experiment. 
In spite of the fact that uo is a steady-state for any choice of e € (0, i), 
numerical dissipation is stronger for larger e as the relative strength of the 
non-linearity decreases. The initial value and the solution to e = 0.1 are 
indistinguishable in the plot. The differ by 0.2% in the supremum norm. 
Continuous piecewise linear functions are not steady-states for the original 
Perona-Malik equation since the non-linearity is not well-defined for such 
functions. In spite of this its numerical counterpart delivers results compa- 
rable to those for small positive e. Figure [4] shows clearly one of the claimed 
enhancements of ( |2.4| over the original Perona-Malik equations. The initial 
condition (in red) 



u (x) 



|10x-5| + 0.2sin(647rx), x € (0,1) 



is evolved to time t = 2 (in blue) with e = 0.3. The new equation can man- 
ifestly differentiate between high low contrast gradients and high contrast 
gradients, which are remarkably well preserved. Observe that the initial os- 
cillatory condition would be left virtually unchanged by the original Perona- 
Malik equation for the same and longer time ranges. Even though this paper 
is concerned with a one-dimensional modification of Perona-Malik, equation 



(2.4) can be considered in a natural two dimensional setting. The qualities 



of the its one dimensional counterpart considered here do carry over to that 
case. Figure [5] shows the evolution of a noisy test image every pixel of which 
has been corrupted by about 15% noise in the gray-scale. For other tests 
and details of the two-dimensional implementation we refer to [9 . 
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Fig. 4. The effect of (2.4 1 on high frequency low contrast oscillations. The initial 



condition (in red) is rapidly evolved to the almost piecewise affine function (in 
blue). 




Fig. 5. The denoising effect obtained for the 2D version of ( |2.4[ | with e = 0.6 and 
Neumann boundary conditions. The initial condition is Lenna's image corrupted 
with about 15% salt and pepper noise. 



6. Conclusions 



A modified Perona-Malik equations has been proposed by strength re- 
duction in the non-linearity via a non-local term. Global well-posedness of 
the new equation and the fact that piecewise constant functions are steady- 
states for it are its main analytical advantages. From the practical point of 
view, the new equation delivers enhanced benefits as compared to Perona- 
Malik and suppresses known shortcomings associated to it. It provides an 
easy, effective, and stable tool for image de-noising. These claims are cor- 
roborated by numerical experiments. 
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